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ABSTRACT 


This  study  examined  the  relationships  between  convective  heat  transfer  and 
barometric  pressure  (P^),  specifically,  how  hypobaric  pressure  affects  the  convective 
heat  transfer  coefficient  (hc).  yThe  U .  S.  Army  Research  Institite  of  Environmental 

ff\  1  f  /  ('  c  „L  l  r‘ 

Medicine  -  High  Elevation  Simulation  Facility  was  used  to  simulate  ■  five 
environmental  conditions:  at  elevation  of  0  (sea  level),  1520  m  (5,000  ft),  3050  m 
(10,000  ft),  4570  m  (15,000  ft),  and  6100  m  (20,000  ft).  In  the  chamber,  constant 
temperature  (30°C)  and  wind  speed  (0.5  m/s)  were  maintained.  The  convective 
heat  transfer  coefficient  was  determined  using  the  naphthalene  sublimation 
technique.  Circular  naphthalene  disks  were  affixed  to  six  sites  on  a  stationary 
copper  manikin.  The  amount  of  naphthalene  weight  loss  through  sublimation  was 
translated  to  h£  using  the  Chilton-Colbum  j-factor  analogy  between  heat  and  mass 
transfer.  As  elevation  increases,  decreases;  hc  should  decrease  accordingly, 
pointing  to  a  diminished  convective  heat  transfer  mechanism.  Hitherto,  the 
relationship  between  hc  and  P^  was  thought  to  be  a  nonlinear  power  function.  Our 
results  showed  that  convective  transfer  indeed  decreased  at  higher  elevation,  but 
hc  and  P^  exhibit  a  linear  relationship.  A  linear  relationship  was  also  obtained 
from  a  theoretical  derivation.  These  results  suggest  that  in  a  hypobaric 
environment,  the  convective  transfer  process  becomes  even  more  diminished,  and 
the  potential  convective  heat  loss  is  smaller  than  the  original  nonlinear  theory 
would  predict. 


Index  terms:  convective  heat  transfer  coefficient,  barometric  pressure,  heat-mass 
transfer  analogy,  Chilton-Colburn  j-factor 


Heat  exchange  between  the  human  body  and  the  environment  not  only 
depends  on  the  more  familiar  factors  such  as  temperature,  wind  speed  and 
humidity,  but  also  is  affected  by  barometric  pressure  (P^).  The  interaction  between 
convective  heat  transfer  coefficient  (hc)  and  has  received  little  attention,  as 
evident  by  the  limited  information  available  in  the  literature.  While  it  is  known 
that  as  elevation  increases,  P^  decreases,  hence  hc  also  decreases  resulting  in  a 
reduced  convective  heat  transfer  process;  the  exact  manner  of  this  dependency  is 
uncertain.  This  study  looked  at,  both  experimentally  and  theoretically,  the 
relationship  between  convective  heat  transfer  and  barometric  pressure,  specifically, 
how  hypobaric  P^  affects  hc. 


METHOD 


The  U.  S.  Army  Research  Institute  of  Environmental  Medicine  (USARIEM) 
High  Elevation  Simulation  Facility  /(altitude  chamber)  was  used  to  simulate  five 
environmental  conditions:  at  sea  level  (760.0  mmHg),  at  elevation  of  1520  m 
(5,000  ft,  632.4  mmHg),  3050  m  (10,000  ft,  522.8  mmHg),  4570  m  (15,000  ft, 
429.1  mmHg),  and  6100  m  (20,000  ft,  349.5  mmHg).  For  all  five  environments, 
the  chamber  ambient  temperature  (Ta)  was  maintain  -  at  30°C,  relative  humidity 
20%  (dew  point  «5°C).  This  high  temperature  was  necessary  to  facilitate  the 
naphthalene  sublimation  process.  At  30°C,  most  human  subjects  will,  to  a  greater 
or  lesser  extend,  perspire.  Once  perspiration  occurs,  it  invariably  contaminates  the 
naphthalene,  and  renders  any  naphthalene  weight  measurement  useless.  To  avoid 
the  problem  of  perspiration  contamination,  a  life-size  (1.68  m2  surface  area) 
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stationary  copper  manikin  was  placed  in  the  altitude  chamber  to  model  the  surface 
area  dimension  and  the  passive  dry  heat  loss  of  a  human  subject.  Since  we  were 
interested  only  in  the  convective  heat  loss,  a  manikin  also  eliminated  any 
evaporative  contribution  generally  involved  with  human  subjects.  A  circular  fan 
aided  in  regulating  the  regional  air  velocity  to  between  0.4  and  0.6  m/s,  depending 
on  different  body  segments  of  the  manikin.  These  levels  of  air  velocity  (see  Table 
1)  must  be  maintained  throughout  all  five  hypobaric  environments. 

Circular  naphthalene  disks  were  attached  to  six  sites:  upper  arm,  lower  arm, 
thigh,  calf,  chest,  and  back,  on  the  manikin.  The  Figure  1  schematic  shows 
locations  of  five  of  the  six  naphthalene  disks.  Placement  of  the  back  disk  mirrored 
exactly  the  chest  disk.  The  naphthalene  disks  (except  chest  and  back  disks)  were 
appropriately  curved  to  conform  to  the  corresponding  body  segment  curvature. 
Airflow  was  directed  normally  at  the  disk  surface.  Since  the  naphthalene  disk 
conformed  to  the  body  segment  curvature  and  sat  directly  over  the  specific  body 
site,  the  local  hc  over  the  specific  site  was  measured,  rather  than  an  average  hc 
for  the  entire  body  segment.  The  regional  air  velocity  and  temperature  were 
measured  at  the  six  disk  sites  using  an  anemometer  tree  system  (TSI  Tree  System, 
see  Figure  2).  The  tree  system  supports  six  pairs  of  omnidirectional  thermal 
anemometers  (TSI  model  1620)  and  thermistor  temperature  probes  (Yellow  Springs 
Instrument  model  705).  Each  anemometer  and  thermistor  pair  was  placed 
approximately  2  cm  above  and  2  cm  away  from  the  naphthalene  disk,  in  such  a 
way  as  not  to  disturb  the  impinging  airflow  to  the  disk.  The  omnidirectional 
anemometers  have  a  range  of  0  -  3  m/s,  and  a  response  time  of  two  seconds. 
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Naphthalene  disks  were  casted  using  scintillation  grade  naphthalene.  With 
a  melting  point  of  approximately  80°C,  the  naphthalene  was  melted  and  then 
poured  into  aluminum  disk  casting  cassettes.  The  naphthalene  quickly  hardened 
at  room  temperature,  after  which  a  nonadhering  cover  plate  was  removed  to  reveal 
a  smooth  naphthalene  surface.  Immediately  after  casting,  the  disks  were  stored 
individually  in  airtight  containers.  All  naphthalene  disks  were  allowed  to 
equilibrate  in  the  chamber  at  30°C  for  24  hours  before  using.  The  casting 
cassettes  were  appropriately  curved  to  conform  to  surface  curvatures  of  the 
upper  arm,  lower  arm,  thigh,  and  calf  of  the  manikin.  Cassettes  used  for  the  chest 
and  back  were  flat.  For  the  study,  the  disks  were  first  mounted  into  vinyl 
retainers,  which  were  in  turn  fastened  to  the  manikin  using  Velcro  straps.  The 
disks  were  weighed  immediately  before  and  after  an  experiment  on  a  balance 
sensitive  to  ±0.01  mg.  The  duration  of  each  experiment  was  55  minutes.  The 
naphthalene  weight  loss,  from  55  minutes  of  sublimation,  were  translated  to  hc 
data  employing  the  Chilton-Colbum  j-factor  analogy  between  heat  and  mass 
transfer  [1], 

Thoery:  heat-mass  transfer  analogy  and  the  effect  of  barometric  pressure 
Heat-mass  transfer  analogy  of  a  sublimation  substance  has  been  traditionally 
used  to  accurately  predict  the  forced  convective  heat  transfer  coefficient  [5,10]. 
Naphthalene  very  conveniently  sublimates  at  room  temperature  and  thus  has  been 
used  by  a  number  of  investigators  to  experimentally  determine  hc  [2,6,11]. 
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Mass  transfer  of  naphthalene  sublimation  hj^  can  be  expressed  as  [6] 


hm  -  R  •  Ta  •  m  /  (Ps  -  Pa> 


U) 


Assuming  the  heat  of  sublimation  is  negligible,  Pg  may  be  considered  as  equal  to 
the  saturated  vapor  pressure  at  Ta  [9], 

log  Pc  =  11.55  -  3765/T 
10  s  a 

or, 

P  _  (11.55-3765/TJ  f9, 

rs  -  10  a  VI 


The  Chilton-Colbum  j-factor  analogy  [1]  can  be  described: 


for  heat  transfer  jV  = - (Pr)^,  with  Pr  = 

11  P*Cp*U  r 

PxM.n’hjjj  p 

for  mass  transfer  j  = - (Sc)  ,  with  Sc  = 


cp’P 


Jm 


u 


P*D. 


(3} 


(4) 


Equating  heat  and  mass  transfer  j-factor, 


p*cp»u 


'C  (Pr)2«  .  P“  n'hm  (Sc)2/3 


15) 


u 


Sc  \2/3 


hc  =  P  *  cp  *  Pam^  *  Nn 


16} 
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hc  =  (  P  • 


(7) 


In  Equation  7,  only  p  and  Dy  are  functions  of  barometric  pressure. 


from  [3  j 


from  [6] 


P 


1.2929 


273.15 


a 


pb 

760.0 


kg/m3 


0.0185  ,  Ta  9 

-  •  ( - ) 

3600  273.15 


760.0 


m2/s 


(8) 

(9) 


Combining  Equations  7  through  9,  and 
properties  at  30°C: 

K  =  0.02636  W*m  1»K1 
cp  =  100 6  J.kg'.K1 

p  =  1.869*  10'5  kg^m'^s1 

Pr  =  — - =  0.7133 

K 

PAM.n  =  0.319 


substituting  in  some  typical  physical 

Dy  =  6.330*  106  m2/s 
p  =  1.165  kg/m3 


Thus, 


,  273.15  r/o 
hr  =  1035.7  •  ( - - - r6 


a 


*  Nn 


760.0 


(10) 
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The  relationship  between  hc  and  P^,  previously  has  been  considered  as  a 
nonlinear  power  function.  Gagge  and  Nishi  [4]  proposed  by  rational  analogy  that 
barometric  pressure  should  change  hc  by  a  factor  of  (P^/760)0'55,  with  P^  expressed 
in  units  of  mmHg.  The  above  theoretical  derivation  of  Equation  10  very  clearly 
shows  that  hc  varies  linearly  with  (Pj^/760),  rather  than  as  a  power  function.  This 
study  also  seeks  to  corroborate  this  linear  relationship  finding  with  experimental 
data. 

RESULTS 

Table  1  gives  the  regional  air  velocity  (±  dispersion  range)  at  the  six 
naphthalene  disk  locations:  upper  arm,  lower  arm,  thigh,  calf,  chest,  and  back,  on 
the  stationary  manikin  at  sea  level  setting.  Figure  3  shows  the  local  air  velocity 
at  the  six  naphthalene  disk  sites,  at  the  five  P^:  760.0,  632.4,  522.8,  429.1,  and 
349.1  mmHg. 

Figure  4  gives  the  local  hc  at  the  six  naphthalene  disk  sites,  at  the  five 
altitudes.  The  hc  was  computed  using  Equation  10.  For  each  experimental  run, 
hc  was  averaged  over  the  55  minute  experiment  period.  The  results  in  Figure  4 
are  the  average  of  19  runs.  A  more  detailed  set  of  the  data  presented  in  Figure  4 
is  tabulated  in  Table  2,  where  hc  is  shown  with  its  corresponding  standard 
deviation  value.  Figure  5  shows  the  same  information  as  Figure  4  in  line  plots, 
to  accentuate  the  linear  characteristics. 

Analysis  of  variance  (ANOVA)  was  performed  on  the  hc  values  at  each  disk 
site,  to  see  if  the  observed  difference  was  indeed  statistically  significant.  The  f- 
values  from  the  ANOVA  tests  are  listed  in  Table  3.  A  /-value  >  2.47  indicates 
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that  the  h£  sequence  does  show  a  significant  difference,  with  p<0.05.  The  Tukey 
test  was  then  used  to  determine  exactly  where  the  differences  lie,  and  the  results 
are  given  in  Table  4. 

The  Tukey  test  was  evaluated  between  hc  at  each  adjacent  elevations.  The 
0-1,  1-3,  3-4,  and  4-6  columns  give  the  test  comparison  between  elevations  of 
0  and  1520m,  1520m  and  3050m,  305pm  and  4570m,  and  4570m  and  6100m, 
respectively.  From  the  results  of  the  Tukey  test  in  Table  4,  for  each  disk  site  the 
difference-between-means  are  greater  than  the  respective  critical-difference- value, 
indicating  that  the  hc  decrements  between  adjacent  elevation  levels  are  all  highly 
significant  (p<0.05). 

DISCUSSION 

The  sea  level  regional  air  velocities  at  the  six  disk  sites,  as  shown  in  Table 
1,  represent  the  base  line  air  velocity  level  which  must  be  maintained  throughout 
all  five  hypobaric  environments.  Unless  constant  air  velocity  is  maintained, 
comparisons  of  hc  from  the  five  environments  become  meaningless.  In  our  study, 
a  circular  fan  was  used  to  aid  in  achieving  this  base  line  air  velocity  level  since 
the  present  altitude  chamber  does  not  stringently  monitor  the  chamber  airflow. 
The  ranges  of  the  regional  air  velocities  are  displayed  in  Figure  4  with  air  velocity 
(in  m/s)  on  the  ordinate,  and  (in  mmHg)  as  the  abscissa.  Data  for  all  six  disk 
locations  are  presented  together  because  all  six  sets  of  data  display  one  single 
result:  constant  air  velocities  were  maintained  throught  the  five  environments.  In 
Figure  4,  at  each  of  the  six  disk  sites,  bars  of  equal  height  signify  that  constant 
regional  wind  speed  was  indeed  maintained  for  all  five  elevations.  Again,  constant 
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regional  air  velocities  ensure  meaningful  comparison  of  convective  heat  transfer 
coefficient  hc  data. 

The  six  h£  bar  graphs  in  Figure  5  were  plotted  with  hc  (in  Wtu^C'1)  on  the 
ordinate,  and  again  (in  mmHg)  as  the  abscissa.  Again,  all  six  graphs  are 
presented  together  to  highlight  one  single  thesis:  h£  decreased  steadily  as  the 
barometric  pressure  decreased  from  760  mmHg  to  349.5  mmHg.  The  decreasing 
trend  of  h£  was  in  fact  a  linear  one.  This  is  accentuated  in  Figure  6  where  instead 
of  bar  graphs,  h£  were  plotted  as  line  graphs.  In  all  cases,  the  decreasing  trend 
of  hc  followed  a  straight  line.  To  further  ascertain  that  the  observed  linear 
relationship  indeed  exists  between  and  hc.  a  regression  analysis  was  performed 
using  the  complete  h£  data  set  of  Table  2.  Table  5  gives  the  results  of  the 
regression  analysis.  For  all  six  disk  sites,  a  linear  equation  best  describes  the  data 
with  correlation  >0.99.  The  conclusion  must  be  that  h£  and  P^  exhibit  a  linear 
relationship. 

The  nonlinear  power  function  dependency  concept  appears  to  be  a  residual 
influence  dating  from  human  calorimetric  studies  of  hc  and  air  velocity,  where  hc 
has  been  shown  to  vary  as  a  power  function  of  air  velocity  [6,7,8].  Most  likely, 
realizing  that  hc  should  vary  with  P^,  and  as  a  reasonable  approximation  the 
power  function  idea  was  adopted.  The  data  in  Figures  4  and  5,  the  theoretical 
derivation  of  h£  in  Equation  10,  and  the  linear  regression  analysis  results  in 
Table  5,  very  clearly  demostrated  this  power  function  concept  to  be  groundless  and 
indeed  the  valid  approach  is  that  hc  and  P^  exhibit  a  linear  dependency. 
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The  implication  of  the  above  results  is  that  at  a  given  hypobaric  environment, 
the  convective  heat  transfer  mechanism  becomes  even  more  restricted,  particularly 
for  human  dry  heat  exchange,  than  was  previously  thought.  As  an  example, 
consider  h£  at  an  elevation  of  4570  m  (P^  =  429.1  mmHg),  and  for  simplicity 
assume  all  other  factors,  e.g.  temperature,  humidity,  remain  constant.  From 


a  power 
function  : 
relationship 


Pb  \Q.55  _  /  429  1  \0.55 
760.0  '  K  760.0  J 


73.0%  of  sea  level  hc  ; 


a  linear  relationship: 


760.0 


-)  =  ( 


429.1 

760.0 


=  56.5%  of  sea  level  hc  . 


The  power  function  overestimation  of  h£  in  this  example  is  16.5%,  i.e.  the 
convective  heat  loss  is  16.5%  less  than  previously  believed.  This  overestimation  of 
hc  becomes  more  pronounced  as  the  altitude  increases. 

When  the  results  of  this  study  are  applied  to  human  clothing  thermal 
resistance  properties,  it  appears  that  at  higher  elevations,  the  intrinsic  insulation 
property  of  a  multilayered  cold  weather  clothing  system  could  well  be 
underestimated  because  convective  heat  loss  is  lower  then  previously  proposed. 
Hence,  for  example,  by  dressing  according  to  the  original  conventional  nonlinear 
power  function  prediction,  an  exercising  person  may  well  be  overdressed  and  suffer 
from  too  much  heat  storage  rather  than  from  the  inherent  cold  of  high  altitude. 
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Equation  10  applies  implicitly  also  to  hyperbaric  conditions  although  no 
specific  hyperbaric  data  are  presented  in  this  report.  At  hyperbaric  P^,  it  can  be 
reasoned  that  h£  would  be  generally  underestimated  using  the  invalid  power 
function  assumption.  Such  underestimation  could  have  rather  significant 
consequences,  since  in  a  hyperbaric  environment,  convection  becomes  the  one 
dominant  heat  transfer  mechanism  [12]. 
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APPENDIX 

Symbols  and  Terminology 

cp  =  specific  heat  (J’kg'^K'1) 

Dv  =  mass  diffusivity  (m2/s) 

h£  =  convective  heat  transfer  coefficient  (W*m'2*K’) 

1^  =  naphthalene  mass  transfer  coefficient  (m/s) 
j'k  =  Chilton-Colburn  analogy  j-factor  for  heat  transfer 
j  =  Chilton-Colburn  analogy  j-factor  for  mass  transfer 
m  =  naphthalene  sublimation  weight  loss  rate  (kg*m  2»s  ’) 

Pa  =  naphthalene  vapor  pressure  in  air  (assumed  =  0) 

=  barometric  pressure  (Torr  or  mmHg,  1  mmHg  =  133.3  Pascal) 

PAM  n=  logarithmic  mean  density  factor  for  naphthalene  sublimating  in  air  (ND) 
Pr  =  Prandtl’s  number  (nondimensional) 

Pc  =  naphthalene  surface  vapor  pressure  (Torr  or  mmHg,  1  mmHg  =  133.3  Pascal) 
R  =  naphthalene  gas  constant  (0.487  mmHg»m3*kg^»K  *) 

Sc  =  Schmidt’s  number  (nondimensional) 

Ta  =  ambient  temperature  (°K) 
u  =  air  velocity  (m) 
k  =  thermal  conductivity  (W'm'^K1) 
p  =  viscosity  (kg^m'^s1) 
p  =  density  (kg/m3) 
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Table  1 


Regional  air  velocity  on  the  stationary  manikin,  at  sea  level 


naphthalene 

disk 

site 

regional 
air  velocity 

(m/s) 

Upper  Arm 

0.538  ±  0.025 

Lower  Arm 

0.576  ±  0.007 

Thigh 

0.690  ±  0.019 

Calf 

0.649  ±  0.043 

Chest 

0.395  ±  0.033 

Back 

0.231  ±  0.013 

air  velocity  data  represent  as  Means  ±  SD 
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Table  2  Regional  h£  at  760.0,  632.4,  522.8,  429.1,  and  349.5  mmHg 


760.0  mmHg 

(0  m)  i 

! 

! 

j  632.4  mmHg 
,  (.1520  m) 

522.8  mmHg 

(3050  m) 

429. 1  mmHg 

(4570  m) 

349.5  mmHg 

(6100  m) 

Upper  Arm 

6.79+0.319 

5.27*0.260 

4.08*0.199 

3.19*0.127 

2.46*0.137 

Lower  Arm 

7.37+0.262 

5.88+0.279 

4.61+0.198 

3.66*0.139 

2.87*0.154 

Thigh 

7.45+0.324 

6.05±0.321 

4.85*0.161 

3.94*0.169 

3.06*0.140 

Calf 

7.38±0.394 

6.06±0.428 

4.96*0.186 

4.27*0.290 

3.26+0.117 

Chest 

5.76±0.260 

4.37*0.254 

3.36*0.183 

2.61*0.107 

1.99*0.130 

Back 

3. 64  ±0.283 

2.69*0.243 

2.08*0.142 

1.69*0.127 

1.29*0.092 

hc  data  represent  as  Means  ±  SD;  unit  of  hc  is  W/(m2»°C) 
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Table  3  /-value  from  Analysis  of  Variance  (ANOVA  test) 


Upper  Arm 

Lower  Arm 

Thigh 

Calf 

Chest 

Back 

/ 

1088.62 

1190.62 

956.93 

485.65 

1029.73 

414.36 

A  /  >  2.47  indicates  that  significant  difference  exists  (p<0.05) 
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Table  4  Tukey  test  results 


Difference  Between  Means 

Critical 

Difference 

Value 

1 

i  0-1 

1  -  3 

CO 

1 

■  - 

4  -  6 

Upper  Arm 

i  1.520 

1.195 

0.886 

0.730 

0.195 

Lower  Arm 

i  1.497 

1.263 

0.953 

0.795 

0.194 

Thigh 

1  1.403 

i 

1.202 

0.912 

0.874 

0.209 

Calf 

1.318 

1.101 

0.688 

1.011 

0.271 

Chest 

1.386 

1.008 

0.758 

0.615 

0.174 

Back 

0.947 

0.610 

0.398 

0.392 

0.169 

column  0-1  compares  data  between  0  and  1520m  elevations 

column  1-3  compares  data  between  1520m  and  3050m  elevations 

column  3-4  compares  data  between  3050m  and  4570m  elevations 

column  4-6  compares  data  between  4570m  and  6100m  elevations 


Difference  is  statistically  significant  (p<0.05)  when  the 
DifFerence-Between-Means  is  greater  than  the  Critical-Difference-Value 
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Table  5  Linear  regression  equation  between  hc  and 


disk 

site 

linear  equation 
from 

regression  analysis 

i 

correlation 

(p<0.0005) 

Upper  Arm 

y  =  0.01054  x  -  1.32 

r  =  0.998 

Lower  Arm 

y  =  0.01100  x  -  1.05 

r  =  0.999 

Thigh 

y  =  0.01065  x  -  0.67 

r  =  0.999 

Calf 

y  =  0.00979  x  -  0.09 

r  =  0.998 

Chest 

y  =  0.00914  x  -  1.30 

r  =  0.997 

Back 

i 

y  =  0.00562  x  -  0.75 

r  =  0.994 

In  the  resulting  linear  equation  y  =  mx  +  b: 

y  =  hc  data  in  W/(m2»°C)  (as  dependent  variable) 
x  =  Pu  in  imnHg  (as  independent  variable) 
m  =  slope  of  straight  line 

b  =  constant  (intercept  value  on  the  ordinate  axis) 
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Figute  4  Local  hc  data 

Figure  5  Local  hc  data  as  line  plots 
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